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We argue that the entanglement Chern number proposed recently is invariant under the adiabatic 
deformation of a gapped many-body groundstate into a disentangled/purified one, which implies a 
partition of the Chern number into subsystems (disentangled Chern number). We generalize the 
idea to another topological number, the Z 2 Berry phase for a system with particle-hole symmetry, 
and apply it to a groundstate in a weak topological phase where the Chern number vanishes but 
the groundstate nevertheless has edge states. This entanglement Berry phase is especially useful for 
characterizing random systems with nontrivial edge states. 


A quantum many-body groundstate can be regarded as 
a mixed state if a system is divided into several pieces and 
some of them are traced out. This enables us to define the 
entanglement entropy and spectrum (Hamiltonian)^ - — . 
These are widely accepted as new tools to characterize 
quantum many-body states. Recently, the entanglement 
spectrum has been successfully applied to topological in¬ 
sulators through the study of edge states along the fic¬ 
titious boundaries between a partition of a system— - — 
This may be reflected by the surprising universality of 
the bulk-edge correspondence— 

Instead of a partition with definite boundaries, an ex¬ 
tensive partition has been introduced: it has been argued 
that the entanglement spectrum of a bulk subsystem 
can be gapless, at which a topological phase transition 
occurs i 15 ' 16 Not only the spectrum of the entanglement 
Hamiltonian but also the corresponding eigenstates are 
useful to characterize the topological phases of the origi¬ 
nal model. We have shown that, taking an extensive but 
asymmetric partition, if the entanglement Hamiltonian 
is still gapped, we can define the entanglement Chern 
number— The entanglement spectrum and entropy of 
the bipartition reflect the topological properties of the 
bulk through the edge states. Instead, the entanglement 
Chern number reflects the bulk property directly. More 
recently, it has been reported that a random partition for 
a translationally invariant system describes a disorder- 
driven topological transition— Thus, the various kinds 
of partitioning for a single pure bulk groundstate reveal 
its topological properties in different environments. 

In this paper, we argue that the entanglement Chern 
number and its generalization to other topological num¬ 
bers are invariant under the adiabatic deformation to 
make the subsystems disentangled. In other words, the 
original entangled state is eventually modified to the sin¬ 
gle tensor product of the two subsystems. In this sense, 
the entanglement topological numbers may be called topo¬ 
logical numbers for a disentangled groundstate or simply 
disentangled topological numbers attached to subsystems. 
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This also has the meaning of a partition of the topologi¬ 
cal numbers. The entanglement entropy and spectrum of 
the bipartition reflect how the states in subsystems are 
entangled in the groundstate wavefunction. On the other 
hand, the entanglement topological numbers clarify the 
property that remains if the entanglement is eliminated 
or disentangled. This process of disentanglement may be 
considered as a purification of the mixed state. After de¬ 
scribing the general idea and the validity of the entangle¬ 
ment Chern number, we discuss the entanglement Berry 
phases applied to the weak topological (WT) phase ) 19 : 20 
which is topologically nontrivial in spite of a vanishing 
Chern number. 

Schmidt decomposition: Let |G) be a many-body 
groundstate of a fermion system, and let A and A be 
a partition of the total system A + A. Then, |G) can be 
Schmidt-decomposed into 

\G)=J2 D s*\®s)®\^), ( 1 ) 

S,X 

where |$ s ) and |5> x ) are, respectively, orthonormal basis 
states for A and A. The normalization of |G) requires 
(G|G) = DsxDtx — tr DD^ = 1. The singular value 
decomposition for D, D = UKV\ where U s i and V x ( are 
unitary matrices and A = diag(Ai, A 2 , • • • , A m , 0, • • • ,0) 
with \t, > 0, leads to 

|G) = X]A/I$<)®l^), ( 2 ) 

t 

where |$<) = \® s )U s t and |$*) = J2 X \^x)V* t . The 
number m of nonzero eigenvalues is called the Schmidt 

number. The normalization condition for |G) is = 

1 . 

Reduced density matrix: The density matrix of the pure 
state |G) is p tot = |G)(G|. Tracing out A and A, we have 
the reduced density matrix p = tr^ p tot in subsystem A 
and its complementary density matrix p = tr^i p tot in 
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subsystem A such that 

P = E = E IWK'M, 

s,t i 

p=Y j i^x^^vi = E ( 3 ) 

x,y L 

The same Af appear in these equations because D^D 
and DD t have the same eigenvalues except for the zero 
eigenvalue. 

Non-interacting fermion system: Let H be a Hamilto¬ 
nian defined on a lattice by 

H = J24hlfc j , ( 4 ) 

where i,j denote some internal degrees of freedom as 
well as the sites. Let 'S’jn be the jth component of the 
nth eigenstate of the Hamiltonian h tot , Yjh t £f t 'S>j n = 
where E is the energy eigenvalues E = 
diag(ei, e 2 , • • • )• The orthogonality and completeness of 
the eigenstates are expressed by Y :l = b mn , 

Y n = Sij, or simply = 'FT' 1 ' = 1, if T is 

regarded as a matrix. Let us define the normal mode op¬ 
erator d n = = Yj c j^j n - Then, the Hamil¬ 

tonian is diagonal, H = Y v end^dn, and the groundstate 
is given by |G) = JX n<np djjO), where uf is a state index 
below which all the states are occupied. 

Let us discuss p and p in Eq. ([3]) for a non-interacting 
fermion system. To this end, assume that all the sites 
and/or internal degrees of freedom are divided into two 
subsystems A and A, which have Na and Na dimensions, 
respectively. They are denoted by a, b € A and a,b £ A 
with a, b = 1,2, • • • Na and a, b = 1, 2, • • • , Na- Define 

P = e~ n /Z, p = e~ n /Z , (5) 

where Z = tr a e -W and Z = tr A e ~ U - For the time 
being, we restrict our discussion to p. Since the entan¬ 
glement Hamiltonian for a non-interacting fermion sys¬ 
tem is also non-interactingjSL we set T~L = Yab c ababCb, 
where a,b € A. Let us diagonalize the Hamilto¬ 
nian h as J2b h °‘bi’bn = where £ = 

diag(ei,£ 2 , • • ■ ,ejv^). Then, introducing the normal 
mode operator f n = Ya{^)na.c a = Ya^Zn, we have 
Td = Yn £ nfnfn- p is now written as 

— £n/ t fn ^ 

P = T-r | _ e \ = [|lra)£n(ln| + |0 n )(l — £ra)(O n | , 
lln(l + e e Y J 

( 6 ) 

where |l ra ) and |0„) are, respectively, the occupied and 
vacant states of the nth fermion defined by /„|0 n ) = 
0 and I 1 ™) = /l 10„), and is the Fermi distribution 
function 


f, as well as £, is often called the entanglement spectrum 
for convenience. To rewrite p in Eq. © in the form 
of Eq. ©, let us define the many-fermion state in the 
occupation number representation \l) = \£\ - ■ • i n - • • (-n a ) 
with the occupation number of the nth fermion l n = 0,1 
and 

Na 

a? = na-^) 1 -^". (s) 

n— 1 

Then, p in Eq. © is now expressed as p = Ye |-0A 2 (tj. 

Correlation matrix: The two-point correlation ma¬ 
trix is useful as an alternative to the entangle¬ 
ment HamiltonianiSi Noting the relation c\cj = 
Y n m '^in^n'^’jmdm, we see that the one-particle correla¬ 
tion function is given by 

C ij = {G\c\c j \G) = P ji: (9) 

where P jt = Y n <n F = Y n <n F is the 

projection operator to the groundstate. We can now de¬ 
fine the correlation matrices in subsystems A and A as 
follows by simply restricting the sites and/or internal de¬ 
grees of freedom in A or A: 

Cab = Cab — Pba , Cgj, = Cab — Pba- (10) 

Alternatively, noting the relationship c' a Cb = 
Yn,m We obtain 

C ah = {G\c\c b \G)=tr\G){G\c\c h 

= (teAPflfm) = (V^Ea’ (11) 

n,m 

where 2 is the diagonal matrix 2 = diag(£i,£ 2 , • • • , £,n a )- 
Thus, the f are the eigenvalues of C. 21 

The complementary reduced density matrix in a 
fermionic representation and the correlation matrix are 
calculated similarly. Solving the one-particle eigen¬ 
value equation YbKbi’bn = Ym^m^mn, where £ = 
diag(£i,-- - ,£,n a ), the entanglement Hamiltonian % = 

Ys. b c a hab c b can t> e expressed in terms of the normal 
mode operator f n (n = 1, ■ • • ,N A ) as U = Y„ £ nflfn- 
Then, p is written in terms of \t) and the corresponding 
A| as p = Ye K)A|(ij, where |£) is the occupation num¬ 
ber representation of f n fermions and A 2 is similar to Eq. 
© but with the Fermi distribution function f n for the 
energy e n . The correlation matrix is also expressed as 

Cab = (V’ s V’ t )fea> wh ere 2 = diag(£i, • • • , £jvJ- 
It should be noted that, as discussed below Eq. 
and p have the same eigenvalues A 2 , and hence the sets of 
eigenvalues {A^} and {A/} are exactly the same, although 
the dimensions Na and Na are generically different. Let 
us suppose Na < Na for simplicity. Then, this is possible 
only when 2 has the same eigenvalues as 2 except for 0 
or 1. Namely, in an appropriate order, we have 

1 —' (£l, £2, * , ^Na ) j 

S = (6>?2, • • • ,€n a >€n a + 1) • • • >6v.a), ( 12 ) 
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where the extra eigenvalues £jv a +i,--- , £,n a are re¬ 
stricted to 0 and 1. We conclude that C for the larger 
subsystem A has the same eigenvalues as C for the smaller 
subsystem A plus extra trivial eigenvalues of 0 or 1. 

1 




0 

|G)=2 A/ |<E>/)®|<!>/) |G>= |<5j)® lOj) 

FIG. 1: Schematic illustration of a disentanglement deforma¬ 
tion. The spectra of £ and £ are identical except for at 0 and 
1 (e = oo and —oo, respectively), and hence the larger system 
inevitably has extra eigenvalues of 0 and 1. If subsystem A 
includes generic eigenvalues 0 < £ < 1, (i) the other one also 
includes the same generic eigenvalues, (ii) each reduced den¬ 
sity matrix satisfies the grand canonical ensemble with a finite 
weight given by Eq. 0 , and therefore, (iii) the groundstate 
|G) is entangled in the sense that it is composed of multiple 
tensor products of the wavefunctions of subsystems A and 
A. If we can deform the spectrum in the left panel into the 
spectrum in the right one (as if we could take the “zero tem¬ 
perature limit”), |G) can be a single tensor product. This 
process can be considered as a purification of the mixed state 
to the pure state, which results in a disentanglement of the 
groundstate wavefunction. 



Suppose that the spectrum of the entanglement Hamil¬ 
tonian £ has a gap at £ = 1/2, as illustrated in Fig. [Q 
and that we calculate some topological numbers of the 
upper bands. In a generic spectrum, the groundstate is a 
linear combination of the tensor product in the Schmidt 
decomposition as in Eqs. 0 and ©. Then, suppose 
that we deform the spectrum adiabatically, making the 
gap larger, and that we eventually reach an extreme spec¬ 
trum with A = 1, i.e., all states have £ = 1 or 0. In this 
case, the largest eigenvalue, Ai, in the singular value de¬ 
composition {A^} becomes Ai = 1 and the other eigenval¬ 
ues are 0, implying that the groundstate is a single tensor 
product. Therefore, this adiabatic process can be consid¬ 
ered as a disentanglement deformation of the groundstate 
wave function between subsystems A and A. From the 
point of view of the reduced density matrix, the process 
is considered as a purification of the mixed state to the 
pure state. On the other hand, topological numbers cal¬ 
culated using the eigenstates of C and C are expected to 
be invariant in this process since the gap between the up¬ 
per and lower bands never closes. Therefore, such topo¬ 
logical numbers, referred to as entanglement topological 
numbers , reveal the topological properties that are invari¬ 
ant even if the entanglement between A and A is elim¬ 
inated. In this sense, they may alternatively be called 
topological numbers of a disentangled groundstate or sim¬ 
ply disentangled topological numbers. If the groundstate 


can be represented by a single tensor product such that 
|G) = |d>i) (g) |#i), a topological number, such as the first 
Chern number or the Berry phase of |G), is the sum of 
the topological numbers of |d>i) and |l>i). This is indeed 
possible since Ai = 1 in Eq. ©• Otherwise, for generic 
nonintegral A^, it may be difficult to define integral topo¬ 
logical numbers simultaneously for |G), |4>^), and |3>^). 
This also implies that a set of entanglement topological 
numbers for A and A may be referred to as a partition of 
a topological number , provided that the bulk gap of h tot 
remains open in the disentanglement deformation. This 
can be checked by the natural sum rule that the topo¬ 
logical number of the groundstate is the sum of the two 
entanglement topological numbers. Note that assuming 
a finite gap for the entanglement Hamiltonian is in con¬ 
trast to the case with edge states for a bipartition, where 
gapless modes of the entanglement Hamiltonian mainly 
contribute to the entanglement entropy. 

Translationally invariant system: So far, we have used 
the subscripts i,j for some internal degrees of free¬ 
dom as well as the sites. We next consider a system 
with translational invariance. To this end, we replace 
i,j —i ia,j/3, where i,j and a, /? denote the sites and 
species, respectively. On the N d lattice in d dimensions 
with the periodic boundary condition, the fermion op¬ 
erator is now denoted by c a {j) and its Fourier trans¬ 
formation is c a (j) = -^J2 k e - lk ' : ’ c a(k), where V = N d 
and = 2tt/N x integer. For a translationally invari¬ 
ant system, the Hamiltonian given by Eq. © becomes 
KT Kajp = — j ) and its Fourier transfor¬ 

mation is given by h tot (i — j) = y e lk '^~Ah tot (k). 
Then, the total Hamiltonian is separated into k sec¬ 
tors, H = SfeSa/3 c a( fc )^a°/3( fc ) c /3(fc)- The Schrodinger 
equation for a given k is given by J2p h t °p(k)'&p n (k) = 
'I , am(k)E mn (k). We assume that the groundstate is 
insulating and that the fermions are occupied up to the 
n F th band, |G) = U n <n F Tl/ d|,(fc)|0), where the normal 
mode operators are defined by d n (k) = ]T) a c a (fc)’F* ra (fc). 
The correlation matrix in Eq. m is then 

C a p{j,j') = ^J2 eiHj '~ j)p ^k), (13) 

k 

where Pp a (k) = En<n F ^0n(k)^* an (k) is the projection 
operator to the groundstate at a fixed k. 

Example 1; Entanglement Chern number. A typical ex¬ 
ample of the entanglement Chern number is the entangle¬ 
ment spin Chern number— for the Kane-Mele modeb 22 
The Hamiltonian h tot (k) is given by a 4 x 4 matrix due 
to the spin and the bipartite lattice. The Rashba term 
mixes the spins, so that it is basically impossible to de¬ 
fine the spin Chern number simply in the momentum 
space— However, projecting the 4x4 P a p(k) ma¬ 
trix in Eq. m into each spin sector a =j~, j, such that 
P a p —> P<jP a pPa, where P a stands for the projection to 
spin er, we have successfully computed the set of entan¬ 
glement Chern numbers (cf,cj.), which indeed describes 
the spin Hall phase when (c-|-, cf) = ±(1, —1)— Although 
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spin is not conserved in a topological insulator in general 
and the time-reversal symmetry guarantees the vanish¬ 
ing of the Chern number, disentanglement between the 
spins implies a nontrivial entanglement spin Chern num¬ 
ber, which justifies the existence of nontrivial spin edge 
states characterizing the phase. The advantage of the 
topological characterization here is that the idea of the 
disentanglement/purification of the mixed state to the 
pure state is simply extended to correlated electrons with 
an interaction. 

Example 2; Entanglement Berry phase: The Berry 
phase here means a winding number for a one¬ 
dimensional system. We apply it to a two-dimensional 
system, based on the method in Refs . 25 ' 26 , to study the 
nontrivial edge states in a WT phase.Consider an 
N x N square lattice with the periodic boundary con¬ 
dition. Let A be a subsystem with ua ladders. The 
remaining subsystem is denoted as A and is composed of 
= N — ua ladders. For the partitions shown in Figs, 
da) and 0(b), we set A = X and A = Y, respectively. 
Let us consider the case A = X. Since the translational 



FIG. 2: Partition A and A in the case of ua = 1, (a) for 
A = X and (b) for A = Y. 


invariance in the y direction is broken, we regard j y as 
the species, and thus the correlation matrix is denoted 

by 

Ci v u,j'„p{jx,3x) = ~y Ye*-U- (-'jy a ,jyp{k x ), (14) 

where 

C jya j' y p(k x ) = e ik ^ j y- j ^Pp a (k). (15) 

ky 

When j y and j' y are restricted within 1 < j y ,jy < n x 
(ny), the above correlation matrix is C ( C ). We as¬ 
sume that the eigenvalues £, n {k x ) of Cj ya j^p{k x ) have 
a spectral gap, as shown in Fig. [Q This is possible 
in general for an extremely asymmetric partition with 
X X. Let IpA^kj;') = {'4^jyOc,ni{k x ),1pj y oc,n2{kx),‘‘‘) 
be the set of eigenstates with eigenvalues £, n {k x ) > 1 / 2 . 
Then, the entanglement Berry phase for X is calculated 
by 7 jv = Imlog [Jj fc U x (k x )\ , where the U(l) link vari¬ 
able is defined as U x {k x ) = det ily + (k x )ip + {k x + Sk x ) with 


6k x = -/( a, unit of the discrete momentum. Likewise, 
solving the eigenvalue equation for C and/or choosing 
R = Y, we obtain the other entanglement Berry phases 
lx, 1 y, 7 Y- The sets of {"lx, lx) and {lY,ly) are con¬ 
sidered as a real-space partition of the conventional Berry 
phases 7 x {k y ) and 7 y {k x ),— as discussed below. While 
the latter are already partitioned into each k y and k x 
for the pure model, the former partition is advantageous 
when we study disordered systems. 

In what follows, we consider an anomalous Hall effect 
model of the Wilson-Dirac type^Zr - 1 with an anisotropic 
Wilson term. The pure model is defined by 

h tot ( k ) = ter 1 sin k x + t<72 sin k y 

+ 03 [m — b x ( 1 — cos k x ) — b y ( 1 — cosfc y )]. (16) 

This model has particle-hole symmetry and its ground- 
state is characterized by the Chern number. In case of 
an anisotropic Wilson term, an interesting c = 0 phase 
appears, which has edge states. This phase has been re¬ 
ferred to as the WT phase ^ >22 In Figs. [3fa) andOb), we 
show the spectrum of the system on a cylinder belonging 
to the WT phase, in which edge states can be seen only 
in (a). The question is their stability: Unless there are 
specific reasons, such states are expected to be unstable 
against, for example, disorder or interactions. 

For a model with particle-hole symmetry, the Berry 
phase can serve as a Z 2 topological invariant^ Fur¬ 
thermore, with translational invariance, the conventional 
Berry phase (winding number along k x ) "/ x {k y ) can be 
computed. Then, at a certain k y where the particle-hole 
symmetry is enhanced to chiral symmetry, 7 x is quan¬ 
tized as 0 or 7 r. The Berry phase 7 X = 7r for a periodic 
system is a topological invariant for a zero-energy state 
localized at the end of a finite chain. This state forms, 
in turn, an edge state at the boundary parallel to the y- 
axis. Namely, from the fc^-resolved Berry phase "f x (k y ), 
we can predict the edge states at the boundary parallel 
to the y-axis. The Berry phase 7 x (k y ) corresponding to 
Fig. Oa) becomes 7 r at k y = 0, 7 r. However, if the sys¬ 
tem breaks translational symmetry, 7 x {k y ) is no longer 
defined. This is a part of our motivation for proposing 
the entanglement Berry phase. 

Let us start with a model with translational symme¬ 
try. In Figs. He) and[5Kd), we show the entanglement 
spectra for the symmetric partition ua = ua = N/2. It 
turns out that the edge states in the real space are well 
simulated by the entanglement spectrum in a cylindri¬ 
cal partition. To study the stability of these states in 
(a) and (c), we calculate the entanglement Berry phase 
for a minimum subsystem X with nx = 1 and its com¬ 
plement X with nx = N — 1, whose spectra are given 
in Figs. He) and Of), respectively. These are indeed 
gapped, and the entanglement Berry phase is therefore 
well defined. We obtain ( 7 x, 7 x) = ( 777 r) numerically. 
This is in sharp contrast to the trivial c = 0 state of the 
model with ( 7 x, "lx) = ( 0 , 0 ) and to the c = 1 state with 
{lx,lx) = (tlO) or ( 0 , 7 r). 
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FIG. 3: Various spectra of the model given by Eq. Hi) 
with t = 1, m = 1.5, b x = 1, and b y = 0 belonging to the 
c = 0 WT phase. Insets show schematic illustrations of the 
boundary or partition of the system. The top two figures are 
spectra of h tot on a cylinder with boundaries (a) parallel to 
the y -axis and (b) parallel to the r-axis. The next two are 
the spectra of C (c) for the subsystem Y with ny = N/2 
(N = 20) and (d) for X with nx = N/2. The edge states 
in these Egures match the top two figures. The third two 
are the same spectra as in (d) but with (e) nx = 1 and (f) 
nx = 19. In (f), 18 states are degenerate at £ = 0 and 1. The 
bottom two are those for the disordered model with the same 
parameters as above. Randomness is included in the mass 
and hopping such that nrij = m + 5mj, tj : p. = t + 5tj ^ with 
a random distribution 5nij,5tj t p, £ [—0.3, 0.3]. (g) Spectrum 
of C for y with ny = N/2 (N = 10) as a function of the 
twist angle (j> y and (h) spectrum of C for X with nx = 1 as a 
function of the twist angle <f> x . 


Finally, we study the same model with impurities. We 
can define the Berry phase even in such a model by im¬ 
posing the twisted boundary condition and by using the 
twist angle (^ x ,(^y) instead of the momentum (k Xl k y ). 
We then obtain 7 X = 0 mod 27 t for <f) y = 0. This is 
expected because the edge states that cross the zero en¬ 
ergy at k y = 0 and n in the pure model are no longer 
distinguished by k y , and hence "f x = 0 = n + 7 r mod 27r 
is observed, even if the edge states remain. Figure [3f g) 
shows the entanglement spectrum for Y under the sym¬ 
metric partition ny = N/2. At <p y = 0 we see that even 
with disorder, the four states seem degenerate near the 
zero energy, which may originate from the zero-energy 
edge states of the pure model. Here, the entanglement 
Berry phase for this state plays a crucial role when dis¬ 
cussing the stability of the WT phase. Let us calculate 
the entanglement Berry phase for a minimum subsystem 
X with nx = 1 and its complement X with nx = N —1. 
The spectrum of X is displayed in Fig. Hh). Note that 
it is indeed gapped, and the computed Berry phase is 
(7X)7x) = ( 7 r J 7r ) even with disorder. These entangle¬ 
ment Berry phases imply that if one divides the system 
into two pieces X and X and regard them as two one¬ 
dimensional chains, each subsystem has each edge states. 
The natural sum rule 0 = ") x = yjy + 7 ^ = 7 r + it mod 2 - 7 T 
indeed holds. In other words, a partitioning of the Berry 
phase enables us to observe the Berry phase 7 r. 

To summarize, we have argued that the entanglement 
topological numbers are invariant under the disentangle¬ 
ment of an entangled groundstate and that they are topo¬ 
logical numbers attached to disentangled subsystems. In 
this sense, the entanglement topological numbers serve 
as a partitioning of the topological numbers. We have 
introduced the entanglement Berry phase to show the 
stability of the edge states in the WT phase. 
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